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Abstract 

We calculated the coupling coefficients for non-linear, quasi-local oscillatory modes of thin 
accretion disks. We found that several of them are non-zero. Mode coupling is a necessary 
condition for a resonance, and thus our results may be relevant for the recently discussed 
QPO resonance model. 

1 Introduction 

Several galactic black hole and neutron star sources in the low mass X-ray binaries show quasi 
periodic variability in their observed X-ray fluxes. Some of the quasi periodic oscillations 
(QPOs) have high frequencies, roughly equal to orbital frequencies a few gravitational radii 
(re = 2GM/c 2 ) from the central object (with the mass M). The high f requency QPOs oft en 
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der Klisl (|200Ch ). 



come in pairs (z/rj, v\f) of twin peaks in the Fourier power spectra (see e.g. 
Low frequency QPOs, with frequencies of a few Hz, are also observed. 

There is no general agreem ent on a physical mechanism exciting the twin peak QPOs. 
Abramowicz fc Kluzniakl (|200ll ) suggested that they may be due to a resonance in some ac- 
cretion disk oscillation modes. This suggestion is further supported by other observational 
properties of twin peak QPOs. The most important support comes from: 
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Twin peak QPOs are oscillations (waves). In black holes an d neutron stars, a linear correla- 
tion w as found between the high and low frequency QPOs ( Psaltis. Belloni fc van der Klisl 
(jl999l )). The correlation proves that the QPO phenomenon is due to accretion disk oscil- 
lations, and not to kinematic effects like e.g. Doppler modulation of fluxes from isolated 
hot spots: thus, QPOs are "waves", not "particles". 

Twin peak QPOs origin is connected to strong gravity. The frequencies of twin peak kH z 
QPOs in microquasars scale with mass, v ~ 1/M, (|McClintock. fc Bemillardl (2003|)). 
This suggests a relativistic origin of them, 
frequency scales as v ~ c/2GM/c 2 ~ 1/M . 



Indeed, in a strong gravitational field the 



3. Twin peak QPOs are due to a resonance. In all four microquasars in which the twin peak 
kHz QPOs have bee n observed, vy/iA, = 3/2. The c ommensurable frequencies obviously 
suggest a resonance dAbramowicz fc Kluzniakl (|200lh ). 



The relevant "QPO resonance pair" should be a combination of the normal modes of the 
black hole accretion flow oscillations. All normal modes are known in the two mathemat ically 
idealized situations. The first one is the thin disk linear diskoseismology; see reviews by iKato 
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Braes. Arras fc Fragile! (|2006h . 



(200 1n andlWagonerl (11999 ). The second one is the thick disk slender torus linear diskoseismology, 



Although picking up the right pair still involves guess work, several authors independently 
argued that the epicyclic modes or their combination s, may be pa r ticularly interesting in th e 
context of the high frequency twin peak QPOs — e.g. IStellal dl99flh IPsaltis fc NormaJ tod ) , 
Kluzniak W. & Abramowicz M. Al (I2002T). Eatol d2003h. iTitarchukl tooi). 

Abramowicz et al. ( 2006 ). Blaes. Arras k. Fragile ( 20061 ) and Blaes et al. ( 2007 ) h ave st udied 



properties of oscillation modes in slender tori, while Kato (see iKatol (|2003T ). iKatol (|2004l ) and 
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Katol(|2007h ) studied a nonlinear resonant coupling between disk oscillations and its deformation. 



In this paper we study non-linear three-mode coupling of epicyclic and other modes in the thin 
disk case, using a quasi-local approach. We shall discuss similarities and differences between our 
and Kato's work in the last section of this article. 



2 Local Analysis of Thin Accretion Disk in Equilibrium 

By rescaling each physical quantity describing an accretion disk in equilibrium with its typical 
value Q* we get Q = <3/Q*~l. In a thin accretion disk all quantities describing stationary 
states change quickly in the vertical direction {z) and slowly in the radial direction (r). That is 

¥ « ¥■ (D 

dr dz 

By rescaling r and z by their typical values r* and 2* so that f = r/r* ~ 1 and z = z/z t ~ 1 we 
write 

dQ = Q*dQ 
dr r* dr 
dQ = Q+dQ 
dz z+ dz 

where dQ/df ~ 1 and dQ/dz ~ 1. One introduces the thinness parameter -C 1, 



/dQ\,/dQ\ _ z*/dQ\ ,(dQ\ z* = 
\dr"\dz) r^dfJ'ydzJ" 
Considering only a small annulus of the disk, i.e. doing a local approach, it is convenient to 
introduce new variables x and y defined as 

= p^, j- = y- (5) 

We assume the radial and vertical sizes of the ring to be comparable, that is £ ~ z± or 
(£/ r o) — fi£ with I = 0(1). The part of the disk we are interested in corresponds to the ranges 
— i < x < £ and — 1 < y < 1. Note that the radius ro and the disk half-thickness h{ro) = ho can 
be taken as typical scales (r* = ro, = h(ro) = ho) so that (3 = ho/ro- 

In our calculations we only consider terms of the lowest order of (3. The stationary quantities 
which will turn o ut to be important are t he angular velocity Vt = Six and the speed of sound c s . 

According to iKluzniak Kital (|200d ) 

H dz r 3 ' 

where n comes from the equation of state p oc p 1 = p 1+1 / n . By using the boundary condition 
c s {z = ho) = and by integrating equation © it is easy to show that 

(7) 

A typical scale for the speed of sound is c St = uj z {ro)ho/yj2n so that 

(8) 



'hi 




2nrl uj zo h \J r&^hl 



Observing that uj 2 zq = w^(ro) = fl 2 K (ro) = r 3 equation ([H]) becomes 

(9) 
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3 Perturbation Equations 

Next thing t o do is to perturb the equa tion of motion. This is done by the Lagrangian approach 
also used by Nowak fe Wagoner ( 199 ll ). that is 



( ^ + -Vp + V$ - -F, 
Dt p p 







(10) 



where p is pressure, $ is gravitational potential and F v is viscous forces. The Lagrangian 
perturbation is defined as AQ(x,t) = Q(x + £,t) — Qo(x,t) where £ is called the Lagrangian 
displacement. Considering a polytropic fluid the perturbation becomes 



2<-i 



I, 
p 



-V, 



p( 7 - i)(v • + P v i e + evjV'® = <(o + ant) 



(ii) 



whe re q l i is the metric t e nsor and D /D£ = d/dt + v kj \/k is t h e Lag rangian flow derivative, see 
e.g. Friedman Schutd ( 19781 ) and Lvnden-Bell fc Qstriker (1967). The nonlinear terms are 
collected in a n ($,), and a v (£) is the contribution of the viscous forces which are ignored in the 
calculations. When using cylindrical coordinates (r, (f>, z) equation (jlip becomes 



2tr 



d 2 e 
dt 2 



20, 



d^ 



a 2 



5 2 $ 

dr 2 



a. 5 2 $ 



■(v-€) 



del 2 ^(v-0 dc 2 s dc dc 2 s dn 



dr 



dr 



n 



dr dr 



n 



dz dr 



(12) 



d 2 ^ 
dt 2 

1 

r 



20 



d£ r 
~dt 



15$ 

r 5r 



s 50 5r d(j) dz dip 



dc 2 s di ' 



v 
r 



(13) 



d^ 



+ 



5 2 $ 



d 2 $ 



dt 2 drdz 



r + -5-»e*-(v-0-5f 



5c s 2 o5(V-0 



5z 2 



5r 
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5c 2 5£ r 5c 2 5f 2 



(14) 



where the operators are defined as 



dt dt dx k 



(15) 



The stationary quantities appearing in these equations have been calculated. In addition ex- 
pressions for the gravitational terms are needed. Some short and simple calculations yield 



" 5 2 $ ~ 

dz 2 

"5 2 $ 

drdz 



z=0 



0. 



J z=0 



5 2 $ 

dr 2 

1 5$" 

r 5r 



2 = 



z=0 



L0 r 



3n 2 , 



(16) 
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Equations (j 12 j) - (|14[) are now divided by ft 2 . When collecting the zeroth order terms in /? the 
resulting equations become 

d^ 



dt 



+ 2C 



d 2 r 



dx 



1 C1 -2\ 5 V 



d 2 ^ 



i 



n 



J_(l_p2)^ + j^! 

2n 9y <9y 



0, 

0, 
= 0, 



(17) 
(18) 
(19) 



where t = Q*t. This is the same result as was derived by iNowak &: Wagonerl ()199ll ). At this 
level of approximation d/di = d/di + d/dcf) and the variable has been introduced as 



^ dx dy 



4 Lowest-order Linear Equations 



(20) 



The coefficients in the linear equations (|18|) and (|19|) do not depend explicitly on t, (ft and x. 
Therefore, the solutions of these equations can be found on the form 



e 



-\{u)t+m<f>+kx) 
—[(Qt+mrp+kx) 



9(y), 



(21) 



where k is a scaled component of the radial wave- vector. Periodic boundary conditions k = irn r /£ 
are assumed, with n r being an integer wave-number. The wave vector in physical units is given 
by k = kh^ 1 = irn r /£. Using the expressions in (|21|) equation (|18|) yields 



2 2nyg - (1 - y 2 )g' 



f = ikco 



(22) 



2n(a 2 -0 2 )-k 2 to 2 (l-y 2 )' 

When substituting this result into equation (|19p a linear second-order differential equation for 
g(y) is obtained 

(1 - y 2 )g" - yU(y)g' + V(y)g = (23) 
where the coefficients are given by 

2(n + l)(a 2 - O 2 ) - k 2 uj 2 z {\ - y 2 ) 



U(y) 



2n 



2n(a 2 -Lo 2 )-k 2 uj 2 z (l-y 2 ) 



(24) 



V(y) 



2n(a 2 - u 2 )(a 2 - O 2 ) - a 2 Co 2 k 2 (l - y 2 ) 



Anuj 2 z k 2 y 2 



{a 2 - uj 2 )uj 2 



(25) 



2n(a 2 -Cj 2 )-k 2 uj 2 z {l-y 2 y 

Equation (j23H must be completed by the boundary condition stating that the Lagrangian 
pressure variations Ap vanishes at the surface of the unperturbed disk. Since 

1 



Ap 



1 



(26) 



and the vertical profile of the pressure is given by p oc (1 — y 2 ) n+1 , this condition takes the form 

/ r)\ n+1 

lim y>(y) (l-f) =0, 

7— >±1 V / 

2n(a 2 - uj 2 r )g' - k 2 cv 2 fg 



y-*±l 

<p(y) 



2n{a 2 - u 2 ) - k 2 u 2 {l - y 2 )' 



(27) 
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Other boundary conditions are set in the equatorial plane. Considering either even or odd modes 
the two possible conditions imposed on g are 

g(0) = 0, or </(0) = 0. (28) 

4.1 Isothermal Disk 

It is possible to get the equations for vertically perturbed isothermal disk from (|23p if n — ► oo. 
But this limit cannot be applied directly to equation (|23l) . because it is derived under the 
condition h -C r. Here the disk thickness, ho, is defined as the height in the z-direction where the 
speed of sound vanishes, i.e. c s (z = ho) = 0. The sound speed is constant in vertically isothermal 
disks. This leads to that the thickness goes to infinity due to h oc n 1 / 2 . Therefore the disk 
thickness must be redefined. This is done with a new thickness, h* , assuming it remains small 
and approximately constant when n — ► oo. The two thicknesses are related as h* = h/(2qn) 1 ^ 2 , 
where q is a dimensionless constant. With the new definition of the disk thickness the scaled 
lengths x and y also need to be redefined. The new variables x* and y* are related to the old 
ones as x* = {2qn) l / 2 x and y* = (2qn) 1 ^ 2 y. The boundary conditions are the same as before. 
The radial wave-vector k* = k/(2qn) 1 ^ 2 has to be reintroduced as well. Now equation (j23l) 
becomes 



r 2 



fg_ 

dy* 2 

which in the limit n — > oo gives 



(2qn- 



y*U +V 

/2qn J ay* 



i2qn ) 



9 = 



a -i {v 2 ~ ul)(a 2 - uj 2 z ) - qa 2 uj 2 z k 2 

qg - yg +- — =§hr — 9 = o, 



{a 2 - uj 2 )uj 



(29) 



(30) 



where the stars in y and k were omitted. The meaning of q is apparent from the following 
relation valid for n < oo 



1 



q 2n\h* 



h 



u z h 



(31) 



This c an be used t o rewrite equation (|30|) as (as was done by Okazaki et all ( 19871 ) and by Kato 
in e.g. iKatol (|2004h ) 



qg" -yg' + {n z -l)g = 0, 
(a 2 - to 2 - c 2 k 2 )a 2 



{a 2 - uj 2 )uj 2 



(32) 



Note that the case n z = is excluded from the discussion because it violates equation (j22|) . This 
case corresponds to purely horizontal oscillations. 

The factor q in equation (|32l) can be eliminated by choosing it to be equal to one. Then 
equation (1321) will b ecome a Hermite equation. The boundary conditions for this equation is 
dOkazaki etall l|l987l )) 



lim e^ 2 ^ 



g(y) = o, (33) 

remain valid in the case of the isothermal 



y— »±oo 

which replaces condition (|27|) . The conditions 
disk. 

The solution to equation (|32p that satisfies the boundary conditions exists only when n z 
takes a positive integer value. Then the solutions are given by Hermite polynomials. Now that 
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g is known it is possible to get / as well and then the Lagrangian displacement £(x,y,(j)) is given 
by (time-dependence is omitted) 



f = iK-^^H{n z ,y)e-'< m ^\ (34) 
n z a z — uif: 

= 2K-^^H (n z ,y)e-'< m *+ ks \ (35) 
n z <t z — u>£ 

i z = KH(n x -l,y) e" 1 ^^, (36) 

where K is a normalization constant, and n 2 = 1,2,3. . .. Hence, a linear mode is determined 
by three natural numbers: m, n r and n z . For a particular mode n z — 1 describes the number 
of modes in the vertical direction. When n z is even/odd the corresponding mode is symmet- 
ric/antisymmetric under the reflection about equatorial plane. The modes with (n r ,n z ) = (0, 1) 
and (1,2) are analogical to the vertical epicyclic and breathing mode of the slender torus, re- 
spectively. 



4.2 Purely Vertical Oscillations 

Instead of assuming an isothermal disk, i.e. here 7 / 1, equation (|23p can be made easier to 
solve by putting k = 0. Then only vertical oscillations will survive. This implies that equation 
(|24l and equation (f25l) become 

17(f) = 2(n + 1), (37) 



= Mfi^if). (38) 

This changes equation (|23|) into 

(1 - y V - 2(n + 1W + 2n ^ 2 -^\ = 0. (39) 

The solutions to equation (I39p are Jacobi polynomials on the form 5 = P^ ! l '' n \y), where 
N=0,l,2... The first four Jacobi polynomials take the forms 



Pt n \v) = i, 

P[ n ' n \y) = (n + 1% 
,(n,n),_x _ (n + 3)(n + 2) /4n + 10^ 3 



n ' ' (y) = v g ' [^j^r - 2y) • (40) 

The dispersion relation in this case becomes 

a 2 = (iV + l)^^ + l). (41) 

The Lagrangian displacement for vertical oscillations is given by 

C = 
= 

where the time dependence is omitted. 
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5 Nonlinear Pressure Coupling 



In the calcul ations to obtain th e nonlinear pressure coupling coefficients we closely follow the 
formalism of ISchenk et all (|200ll ). A compact way to write equations (|12|) -(|14| ) is 



g + S §f + T«*) = 0, 



(43) 



The solution to equation (I43p can be assumed to be on the form of normal modes 

£ a (x,t) = e- iu <* t t a (*) (44) 
where £ Q is the corresponding eigenfunction. Inserting the assumption into equation (|43p gives 

- uj a 2 £ a (x) - iu a S£ a (x) + T£ a (x) = 0. (45) 

The linear modes are determined by three integer numbers m, n r and n z and form a complete, 
generally non-orthonormal basis in the Hilbert space Ti of all possible Lagrangian displacements 
with the scalar product defined as 



Jv 



(46) 



Inspecting equations (|12p - (ll4l) we find that the operators S and T are res pectively anti - Herm itian 
and Hermitian with respect to the above scalar product. As shown by Schenk et al ( 200ll ) the 
eigenfunctions of the modes satisfy the pseudo-orthogonality relation 



(w a +w/j)(e a ,^) + (e al i%) = o 



(47) 



Looking at equation (|43l) and now including nonlinear terms in £ but still only zeroth order 
terms in (3 the equation will become 



d 2 £ d£ 



(48) 



where represents effects of nonlinear accelerations. Thanks to the completeness of the 

linear- mode eigenfunctions, the solution of the nonlinear equation (|48jl can be found on the 
form 

£(t, 4>, x,y) =Y^ Ca(t)£ a (<t>, x, y) + cc (49) 

a 

where 'cc' denotes complex conjugate and £ a denotes eigenfunctions of the linear modes, a = 
(m, n r ,n z ). The problem with non-orthogonality of the basis can b e handled by carryin g out this 
expansion in phase-space rather than just in configuration space (pchenk et all (|200ll )). Hence 
the time derivative of the Lagrangian displacement is expanded as 



dt£(t, <f), x, y) = ^2 -iu a c a (t)£ a ((/), x, y) + cc. 



(50) 



Considering three different expressions using equation 



(51) 



(52) 
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(Z a ,iS£) = (iS£ a , cpip) = cp(-iS£ a , £ p ) (53) 
and then summing these gives 

(£«, + + iS£) = c/3((€a. C/j) k + w/j) + (»S£ a , ^}). (54) 

Due to that non-zero terms only exist when a = (3 (as follows from the pseudo-orthogonality 
condition) equation (154)) will become 

(S a ,v a e + ii + iSt)=c a (2u a (t a ,S a ) + (S a ,iSZ a )), (55) 

where (2u) a (£ a , £ a ) + {€ a ,iS€ a )) is defined as & a . 

Looking again at equation (|48|) and performing a scalar product with £ a gives 

§ + 50 + (C, re) = <C o(0>. (56) 

' — J — ' 

Expanding I and II separately will give 

d d dc • 

= -i-^c a b a + -^(^ a ,uj a $) = ~ ib a-^- + (£ a ,iw a €), (57) 

= UJ a {u a £ a + iS£ a ,$,) = L0 a {£ a , UJ a $, + iSg), (58) 

where equations (f5"5"j) and (|45p have been used. Equation (f5"6"|) now becomes 

-ib a ^+ Lu a {£ a , i£ + u, a £ + iS£) = {£ a , a(£)>. (59) 
Rearranging equation (foT?|) gives 

+iw Q C Q = ^-(£ a ,a(£)), ( 60 ) 

which can be seen as a set of ordinary differential equations for many harmonic oscillators. The 
scalar product on the right-hand side of equation (f60|) can be expressed as 

(£a> = X] K «/?7 C /3 C 7 + K a^ c } c l + KaPyCpC* + K^jCpC*, (61) 

0,7 

when considering only quadratic terms of £ in a(£) and where k is defined as 

K «/37 = (£a> a (£/5>£ 7 ))- (62) 

In order to get the coupling coefficients it is necessary to find an explicit expression for th e 
acceleration. It can be done by adopting the formalism used by Kumar Sz Goldreich ( 19891 ). 
starting with expanding the Lagrangian density as 

L = L + Li+L 2 + L 3 + ..., L n oc\£\ n , (63) 
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where L\ is the first order contribution describing the stationary flow while the second order 
part, L2, gives the equation for linear perturbations. The higher order terms gives the nonlinear 
evolution of the linear modes. Expanding L leads to that the Euler-Lagrangian equations can 
be applied as follows 

EL(L) = EL(L ) + EL{L X ) + EL(L 2 ) + EL(L 3 ) + . . . = 0. (64) 

Considering only the first orders of perturbations gives 

EL{L 2 ) = -EL(L 3 ) = pad) (65) 

since both EL(Lq) and EL{L\) are zero. a(£) is as before an acceleration due to the nonlinear 
terms. The equation describing the Lagrangian density is 

L = \p{v l + dtC + v k V k Cf - P^-r - p$(x l + C ), (66) 
2 7 — 1 

where v l , p and p are the velocity, density and pressure of the stationary flow, <3? is the gravita- 
tional potential and J is the Jacobian defined as 



J 1 " 1 = 1 + (1 - 7)V£ + i(l - 7)((1 - 7)(V£) 2 - Y,C J Y^) - 

1(1 _ 7 )(( 7 _ l) 2 (V£) 3 + 3(7 - L)(V£)V^*V fc ^ + 2V^V J £ fe V fc £\ (67) 

To get L3 only terms with £ 3 -dependence contribute. This leads to that only the second term in 
equation (fB6^) contributes to L3. Using equation ([6"T|) and the fact that a(£) = EL(L 3 )/p (from 
equation (I65p ) it is now possible to write down an expression for the coupling coefficients. The 
three-mode coupling coefficient K Q/ g 7 is then given by 

KaPl = \l v ^ 7 ~ 1 ) 2??ar?/3r?7 + 3p ( J ~ + 2 PVa^ - pd^VjVjVjfc*} dV, (68) 

where 7 = 1 + 1/n and square bracket denotes symmetrization. The following definitions has 
been done to simplify the equation 

Va^V-Ca, Va? = Vcfr = hC,)- (69) 

While the first three terms of the integrand in equation ([68]) are responsible for the pressure 
coupling, the last term governs the gravitational coupling. In the thin disk p ~ (3 2 p and the last 
term is therfore 6eia-times smaller than the previous ones. 

5.1 Isothermal Disk 

We have then looked at the case of an isothermal disk. The equation for the coupling coefficients 
(equation (|68j) ) will simplify due to the fact that 7 = 1 for isothermal disks and if only terms 
with the lowest order of (3 are kept. The only surviving term in (|68p is then 

= / p(v,e)(v^)(v4)dy, (70) 

where 

(v.exv^xv^) = (d& + + v k jm ^){d k ^ + ru™)- (7i) 
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Because low power of (5 is desired the only term kept from (|7ip is 

mime^dkQ- (72) 

If symmetry in <f> direction is assumed and the indicies are run over the scaled r and z, see 
equation (0), the result will be 

Hdyomime,) + (wcs^xw}- (73) 

Looking at one term at a time and using the Lagrangian displacements from equations (1341) - 
(f36|) gives 

(^CX^X^) = Ae-^ m ^ + ~ k ^H na H np H rii 

n a n^(G^ - u?)(a% - u>%)(g% - Q?) ' 

(^owc^e;) = B^ m ^ k ^u ni R na ^H nfi _ x 

B = K 3 k a k^G 2 ^ 

2n 7 (a| - u> 2 )(a 2 - u; 2 ) ' 

(^)(^)(<%£ 7 ) = Ce-^^^H^H^Hn^ 
K 3 k a k^G^(n 3 - 1) 
C = 2(^-4) ' (?6) 

iK'%k^G}G 2 {n a - 1) 

2n 7 (a 2 -^ 2 )(cr 2 -w 2 )' { " } 
E = K 3 (n a -l)(n p -l)( ni -l), (78) 



{dyOm r p ){d s Q = Fe-^ m ^ +k ^H na ^H n0 H nj 



n^|-^)(a 2 -u> 2 )' 1 J 



G = K ^ n l~ l \ (80) 
2(<r 2 - w 2 ) 
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Here H' m = mH m _i and H' Q = 0), m a p^ = m a + mp + m 7 and /c Q /3 7 = k a + + fc 7 have 

been used. Adding the individual terms together gives 

+DH na -2H n/3 -xH n ^ + EH na _2H n0 -2H ni _2 + FH na -xH np H n ^ 
+GH na „ 1 H np _ l H nj _ 2 + JH^Hn^H^e-^^+^^dV. (82) 

The only (^-dependence is in e - *" 1 "' 37 ^ = cos(m Q!( g 7 (/>) — isin(m aj g 7 ^)). The integral is over an 
even interval, which implies that ra a ^ = otherwise k is equal to zero. This also applies to the 
x-dependence and gives /c Q/ g 7 = 0. For isothermal disks the pressure is p = c 2 poe~v I 2 . Putting 
all this together gives 

AttIc^qPo r°° _-2 /2 
K ^Pi = — ^3^3 — J e {AH na H n/3 H ni + BH na _iH nf} _iH ni 

+C H ricc -iH n/3 _2H ni -i + DH ria -2H ni3 -iH ni + Si/ na _2i^ n(3 -2-ffn 7 -2 
+FH na _iH n/3 H ni + GH na _iH nfj _iH ni -2 + JH not H n0 ^iH ni -i)dy. (83) 

To simplify further we have looked at different properties of Hermite polynomials. The product 
of two Hermite polynomials can be expressed as 

min(i,j) 

HiHj = i\j\ £ -T7-, ^ v .H l+j - 2r - (84) 

Multiplying equation (|84|) with a third Hermite polynomial along with the exponential function 
that appears in equation (|83[) and then integrating yields 

/ H^Hke-y l 2 dy = -±- £ — — / F i+i _ 2r .ff fe e^ ' 2 dy. (85) 

The orthogonality condition for Hermite polynomials is 

/oo 
tfiflj-e-* /2 dy = i!v / 2^^-. (86) 
-00 

Using this in equation (|85l) gives 

J H iHj H k e^ 2 dy = ^ ~i\j\k\ £ r , (i _ r) , ( ■ _ r A w- (87) 



r=0 



This gives the condition that 

/•oo 

/,,/.. = / 7 2 <l/y 7^ 

JO 



11 



S. Fogelstrom, L. Levin & R. Slapak 



when i + j — 2r = k, that is r = (i + j — k)/2. Using this equation (JH7J) becomes 

/7r i\j\k\ 
2 p±EfcJi ^--^ 2 + fc )! ( - i+ 2 7+fc )' 

with the constraints that i + j > k, i + k>j,j + k>i and (i + J + fc)/2 £ Z. Now equation 
(|83l) can be expressed simpler 

^a(3"f — ^3 3 ( ^In a ,rifj,n 1 ~\~ -^-^n a — l,n^ — l,n 7 "I - C-^n a — l,n^— 2,n 7 — 1 ~\~ -D In a —2,np — l,n 1 

p r 

~^~EI na —2,nf-j—2,nj—2 ~i~ ^na-l^^n-, ~i~ GI na ^\^ n ^ — \^ n ^—2 ~\~ J In a ,np — l,n 7 — 1^ • (90) 

5.2 Purely Vertical Oscillations 

As in the isothermal case we only want terms with low order of (3 in equation (I68p . but here 
7 7^ 1. Since £ z is the only surviving term we get 

K 3 

VaVprh = V[aV M = riafr = ^e- 3im ^(dyP na )(dyP np )(dyP ni ). (91) 

P r 

This gives coupling coefficients 

K 3 



1 f K 6 

fh = o / (p(7-l) 2 + 3p(7-l) + 2p)^e- 3 ^n^ J PnJ(5 g P rifl )(9 s P n7 )dy 

/ e- 3 ^ P7 (7 + l)a-P re J(^P n ,)(^P n7 )dF. (92) 



1 



" 2 /?3r 3 iy 

In this case the pressure is p = po(l — y 2 ) n+1 . For the same reasons as in the isothermal case, 
m a/3'y = 0- Using this the final expression will be 

K 3 f l 

*afr = 7«T^07(7 + 1) / (1 " y 2 ) n+1 (dyP na )(dyP np )(d y P ni )dy. (93) 

6 Summary of the Results 

In the case of isothermal disk, see equation (|83j), the following selection rules yield: 



7ij > 1 where i = a, 0, 7 otherwise k is undefined 
n a ,np,nry > 1 gives the following constraints on Iijk- 

(i + i + /e)/2 G Z 

i + j > k 

i + k> j 

j + k > i 

It can be noted that k in this case is always nonzero. Some examples of coupling coefficients 
follows: 



Mil - ^,3 *\N 
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«ii2 = ^ (2A + G + J) 

47dc2po 
^222 = ^^(8^- 



«223 — —3373—0^ y 2 



^233 = ^§^{^A + 6B + 2C + E + 2G + 8J) 

K 333 = ^f (615 + 36^)^ 

In the case of only vertical oscillations, see equation (|93p . the selection rules becomes more 
compact. The only cases when k is nonzero are when n a , np and ra 7 is all odd, or when two of 
them are even and one is odd. 



7 Discussion and Conclusions 



In this work we examined nonlinear interactions between linear modes of the thin accretion 
disk. Our a pproach is analogous to t hose u sed in problems of nonlinear stella r pulsa ti on (see 
e.g. wor ks oflLvnden-Bell & Ostrikerl (Il967h . [Friedman fc Schutd (fl978T ) . ISchut zl (ll980ah . lSchutz 
(Il980bh . ISchenk et all (120011 ) and many others). The Lagrangian displacement that characterize 
the resulting oscillations is given by equation (fl9l) . the time-dependent coefficients c a (t) are 
solutions of the coupled-anharmonic-oscillators equations (|60|) . Hence, a problem of solving 
nonlinear partial differential equations has been reduced to a problem of finding solutions of 
many ordinary differential equations (Schenk et al 2003). In the rest of the paper we concentrated 
on the effects of quadratic nonlinearities and gave explicit formulae for the three-mode coupling 
coefficients K a f3 T 

Our key assumption is a completeness of the eigenfunctions of the linear normal modes. We 
note that previous analysis of linear oscillations by groups of Kato and Wagoner result in such 
complete set of eigenfunctions. For simplicity, we consider here an idealized situation of an 
isothermal disk. Moreover, another substantial approximation is that the nonlinear interactions 
among modes are studied in a small annuli whose radial extension is of the order of the thickness 
of the disk. Under these simplifications the eigenfunctions of the disk are remarkably simple. 
The next necessary step is to use the WKBJ approximation to describe the radial structure of 
the modes. Preliminary results show that the three modes interact only locally, close to radii, 
where a condition k a ± kp ± /c 7 ~ is met (k^ are wavevectors of the modes). 

The coupling coefficient characterize energy exchange between modes due to nonlinear in- 
teractions. In our analysis the disk is made of a polytropic flow with no magnetic field. Hence 
the only contributions to the coupling coefficient are pressure and gravity forces. As we demon- 
strated in section [5] the contribution of the latter is negligible for a thin accretion disk. 

In section O simple selection rules for the three-mode coupling coefficients are derived. They 
are based on the symmetry properties of the modal eigenfunctions. In principle, a nonlinear 
interaction among modes is possible only when (n a + np + n~) is an even number, where 
denote quantum numbers that characterize vertical structure of the modes. The second necessary 
condition is the "triangle inequality" : each quantum number must not be greater then the sum 
of the remaining two. Hence, nonlinear interactions of two oscillation modes (with quantum 
numbers n a and np) may lead to an excitation of many others modes. 
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There is an important physical difference bet ween the typ e of re s onanc e stu died here an d 
that studied in several recent papers by Kato ( see iKatcl l|200a ). iKatol (120041 ) and|KatcJ (|2007l )) 



We deal with a resonance that is excited by a nonlinear three-mode coupling of disk oscillatory 
modes, and Kato consideres a nonlinear resonant coupling between disk oscillations and its 
deformation. In Kato's picture the deformation (assumed to be a fixed perturbation on the 
disks) is essential, as it provides the energy for the resonance excitation. Contrary to what is 
studied in the present paper, in Kato's resonance no direct couplings between the disk oscillations 
are considered. 

The mathematical formalism used in this paper may also be used to study several other issues: 
linear and nonlinear viscous damping of the modes, possible excitation mechanisms, or low- 
frequency modulations of the hig h-frequency oscillations proposed by A bramowicz & Kluzniak 
as an underlying mechanism for iPsaltis. Belloni fc van der Klid (| 19991 ) correlation. They will 



be addressed in our future research. We believe that the formalism may be useful for a future 
development of the general theory of nonlinear diskoseismology. 

Research reported here was done in connection with our master diploma work at the Physics 
Department of Goteborg University, under supervision of M.A. Abramowicz. We thank him for 
all his support and encouragement. We also thank J. Horak for introducing us to the mathe- 
matics of non-linear coupling and for constant help and guidance. We thank P. Rebusco aswell 
for explaining us the multiple scales method. Part of this research was carried out at Nordita 
in Copenhagen and at the Institute of Astronomy in Prague. We thank these institutions, In- 
ternational Astronomical Union, Copernicus Center in Warsaw and also Goteborg University, 
for their support through several research and travel grants, in particular the Polish KBN grant 
N203 009 31/1466. 
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